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In analogy to Glauber's analysis of optical coherence, we adopt an operational approach to in- 
troduce different classes of atomic coherence associated with different types of measurements. For 
the sake of concreteness we consider specifically fluorescence, nonresonant imaging and ionization. 
We introduce definitions of coherence appropriate to them, which we call electronic, density and 
field coherence, respectively. We illustrate these concepts in various descriptions of Bose-Einstein 
condensation, showing that each of these descriptions makes different implicit assumptions on the 
coherence of the system. We also study the impact of elastic collision on the field and density coher- 
ence properties of atom lasers. 

PACS numbers: 03.75.-b, 42.50.Ar, 03.75.Fi, 42.50.Ct 



I. INTRODUCTION 



Quantum optical coherence theory is based on the 
factorization properties of normally-ordered correlation 
functions of the electric field operator |j| . This is a direct 
consequence of the fact that most optical experiments de- 
tect light by absorption, i.e. by "removing" photons from 
the light field. But the situation is not so simple in the 
case of matter-wave fields, and in particular for atomic 
de Broglie waves. This is because atomic detectors can 
work in a number of different ways: For instance, one can 
choose to measure electronic properties of the atoms, or 
center-of-mass properties, or both. Measurements may 
or may not remove atoms from the field, hence the role 
of the annihilation operator is not as central as for light 
fields. 

Due to the added complexity of that situation as com- 
pared to the optical case, no unified theory of atomic co- 
herence exists to-date. The major objective of this paper 
is to present a step toward this goal, extending the ideas 
of optical coherence theory to introduce several types of 
matter- waves coherence, in particular "field coherence" 
and "density coherence." 

In analogy to Glauber's analysis of optical coherence 
|EJ, we adopt an operational approach where different 
classes of atomic coherence are associated with different 
types of measurements. For the sake of concreteness we 
consider specifically three classes of measurements: fluo- 
rescence, nonresonant imaging and ionization. Section II 
briefly reviews the outcome of these measurements, and 
introduces definitions of coherence appropriate to them. 
We call them electronic, density and field coherence, re- 
spectively. In the case of bosonic atoms, field- coherent 
states are easily seen to correspond to the usual Glauber 
coherent states [0. In the single-mode case, density- 
coherent states are simply given by Fock states, but a 
more general discussion is required for the multimode 
case. Section III illustrates these ideas in various de- 
scriptions of Bose-Einstein condensates. Section IV fur- 
ther develops the concept of multimode density coher- 
ence, which is applied to the case of the atom laser in 



Sec. V. In this latter example, we illustrate in particu- 
lar how elastic collisions, while strongly modifying the 
field coherence of the device, have almost no effect on 
its density coherence. Finally, Sec. VI is a summary and 
outlook. 



II. FIELD AND DENSITY COHERENCE 

A physically appealing operational way to discuss 
matter-waves coherence relies on the analysis of specific 
detection schemes, in complete analogy with the optical 
case. As we shall see, this approach naturally leads to 
the need to associate different classes of coherence with 
different types of measurements. In addition, it builds 
useful bridges between concepts familiar in quantum op- 
tics and methods of traditional manybody theory. 

In this paper we consider for concreteness a bosonic 
atomic Schrodinger field described by a multicompo- 
nent field operator ^(r, t), with components $i(r, t),i = 
1,2,... We normally think of the atoms comprising this 
field as adequately described in the Born-Oppcnhcimcr 
approximation, so that r describes their center-of-mass 
motion and the index i labels various electronic states. 
For bosonic atoms, one has then 
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Our goal is to characterize the statistical properties 
of one or more components of this field. One familiar 
method involves the detection of light fields interacting 
with the atomic sample, in the hope that properties of 
the Schrodinger field can be inferred from it. All of laser 
spectroscopy relies on this approach, although it is not 
normally cast in terms of matter- wave fields. Another 
approach, which we will find useful in some respects, 
involves ionizing atoms from the sample and studying 
the properties of the emitted ions or electrons. But this 
method typically also relies on the interaction of the 
atoms with a light field. Hence, a generic Hamiltonian 
describing a measurement scheme for the properties of 
the atomic Schrodinger field is of the form 



1 



V = 



ij 



V$j(r)[d ir E(r,i)]f,(r), 



(2) 



where E(r, t) is the electric field operator, dy is the 
dipole matrix element between electronic states i and 
j, and we have assumed for simplicity that the electric 
dipole approximation gives an adequate description of 
the atom-field interaction. 

We consider the situation where the electromagnetic 
field consists of a classically populated field mode of am- 
plitude Eg, wave vector ko and polarization e, and a se- 
ries of weakly excited sidemodes of wave vectors k^ and 
polarizations eg. In that case, the Hamiltonian (^|) be- 
comes 



V 



E 



O i(ko-r— uot) 



where we have introduced the Rabi frequencies 
ttij,o( r ) = dijE Q (r){eij ■ e)/% 



(3) 



(4) 



corresponding to the classical driving field and the vac- 
uum Rabi frequencies 



Slij,e{r) = di 3 £i{r)(e i: j ■ ee)/h. 



(5) 



In these expressions, is the direction of the Atomic 
dipole of magnitude dij for the i <-> j transition, and 
Si = [HuJe/ZeeV} 1 / 2 is the "electric field per photon" in 
mode I. 



A. Resonance fluorescence 

In resonance fluorescence measurements, one proceeds 
by shining a laser quasiresonant with an electronic tran- 
sition g <-> e, and measuring, e.g., the fluorescence spec- 
trum 



S(u) 



dre 



-(E-(t ,0)E + (t ,t))+c.c. (6) 



at the location ro of a photodetector. In that expression, 
E + (r , t) and E~ (r , t) are the familiar positive and neg- 
ative frequency parts of the electric field operator, and we 
have assumed stationarity to identify the Fourier trans- 
form of the first-order correlation function of the field 
with its spectrum. It is well known that except for unim- 
portant retardation effects, one has that 



E+{r ,t) cx *t( r0jt )$ ( r „ jt ) j 



(7) 



so that 



S(u)<x. / dre-' lUT 

(*t (r , 0)* e (r , 0)*t (r , r)* 3 (r , r)) + c.c. (8) 



Resonance fluorescence has been studied in consider- 
able detail, both experimentally and theoretically. The 
resonance fluorescence spectrum of a two-state atom is 
known to consist of a coherent and an incoherent contri- 
bution. For a perfectly monochromatic excitation laser, 
the coherent spectrum S co h consists of a delta-function at 
the laser frequency, while the incoherent spectrum is the 
famous Mollow three-peak spectrum. Mathematically, 



5 co/l cx / dre™ 

(£t ( ro , 0)# e (r , 0)><#i(r o , r)# 9 (r , r)> +c.c. (9) 

and Si nc (u>) = S(u>) — S co h(^>)- In physical terms, this 
means that coherent effects are associated with the fac- 
torized correlation function 

(*t (r , 0)# e (r , 0)) (*t ( r0j r)* 3 (r , r)>. (10) 

These considerations justify defining the "electronic co- 
herence" of the matter field ^(r, t) in terms of the factor- 
ization properties of normally ordered correlation func- 
tions of the field polarization operator 

E_(r,t) = *t( r)t )$ e ( rjt ). (ii) 



B. Off- resonant imaging 

In contrast to resonance fluorescence, off-resonant 
imaging involves a strongly detuned electromagnetic field 
interacting with the atoms in the sample in such a way 
that it induces only virtual transitions. After adiabatic 
elimination of the upper electronic state of the atomic 
transition under consideration, the interaction between 
the Schrodinger field and the radiation field is described 
to lowest order in the side-modes by the effective Hamil- 
tonian 



V 



d A r 



°o( r )^ atC -i(ko-k,y r 



So 



(12) 



where the atom-field detuning Sq = uj a — ujq is such that 
| Jo | ^> |fio(r)|) an d we have omitted the index label- 
ing the ground state component of the Schrodinger field, 
which can now be considered as scalar. 

There are a number of ways in which off-resonant imag- 
ing can be applied to the determination of specific proper- 
ties of the Schrodinger field. For instance, one can detect 
interferences between the classical incident field and scat- 
tered light, as in the MIT experiments §. This results 
in a signal proportional to the density (p(r, t)}, where we 
have introduced the field density operator 



(13) 
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whose expectation value is the local density of the sam- 
ple. 

Alternatively, one can measure the spectrum of the 
scattered light in a fashion familiar from resonance flu- 
orescence experiments ||. For side-modes initially in a 
vacuum state, the most important nontrivial contribu- 
tion to the fluorescence signal T is proportional to the 
intensity |f2o| 2 of the incident field, 

T= W_y 2 r Mr , f t+At dTdT , 

$0 £ J Jt 

Xe i((k -k 4 ).(r- r ')-( W o-^)(T-r'))^( r)r ) / g( I .' )T /^ ) (14) 

and hence is sensitive to the second-order correlation 
function of the Schrodinger field density. Indeed, it can 
be shown that any measurement involving the electro- 
magnetic field scattered by the atomic sample under con- 
ditions of off-resonant imaging is determined by correla- 
tion functions of p(r,t), 



D (n) ( Xl , ...,x n ) = (p(xi) . ..p(x n )), 



(15) 



where Xi = (r<,ii). 

In analogy with the optical case, we therefore define 
a Schrodinger field as being density- coherent to order M 
if its density correlation functions D^ n \xi : ...,x n ) fac- 
torize for all n < Af. ^From this definition, it is obvious 
that single-mode density coherent states are the familiar 
number states. But the situation is more complex for 
multimode fields, to which we return in Sec. IV. 



C. Ionization 



extracts information about the state of the Schrodinger 
field ^ g (r,t) by standard atomic physics methods, such 
as, e.g., the detection of the quasi- free electrons of the 
continuum states. It is therefore sufficient to describe 
the light field via its (possibly time-dependent) classical 
Rabi frequencies f2j between the levels \g) and \i), so that 
the atom-field interaction reduces to 



V = h 



tfVn^r, i)fl(r)$ 9 (r)e- ,uti + H.c. (16) 



For ground-state atoms cooled well below the recoil tem- 
perature and tightly focused laser beams, the spatial size 
of the atomic wave function is much larger than the laser 
spot and we can approximate the electric field E(r) by 
E(r) ~ E(5(r — r ), so that Eq. (16) becomes 



V = ^^r ,t)#t( ro )* g ( ro ) e — Lt + H.c. (17) 

i 

We assume for simplicity that the center-of-mass wave 
functions of the continuum states of the atoms are well 
described by plane waves of momentum q, so that the 
single-atom Hamiltonian TCq may be expressed as 



= Tin 



where 



Hi, 



tU0i„b] J) 



Here we expanded \&i(r) in plane waves as 



(18) 



(19) 



(20) 



1. Physical model 

The reason resonance fluorescence and off-resonant 
imaging yield signals proportional to correlation func- 
tions of bilinear products of components of the 
Schrodinger field is of course that the electric dipole inter- 
action is itself bilinear in the Schrodinger field operators. 
This raises the question as to whether it is possible to 
measure correlation functions of *P(r, t) itself, as in the 
case of optical fields. This can be achieved if instead 
of making measurements on the radiation field, one de- 
tects the atoms directly One possible scheme that 
achieves this goal is the ionization method that we now 
discuss. 

Consider a detector consisting of a tightly focussed 
laser beam that can ionize atoms by inducing transitions 
from their ground electronic level \g) to a continuum 
level We are interested in measuring properties of 
the ground state component ^ g (r) of this field, which is 
electric dipole-coupled to continuum states ^(r). 

In contrast to the preceding measurement schemes, we 
are no longer interested in the dynamics of the light field, 
whose role is merely to ionize the atoms. Rather, one 



with [bi t q, = S m >5 it ', and uj iq = hq 2 /2M + uji. 

In the following we take the atomic system to be ini- 
tially in the state 



= KV'i.q},^)- 



(21) 



To first order in perturbation theory and for CW beams, 
the transition probability away from that state during 
the time interval At is 



w ~ I^MI 2 / dr / dr' 

x({^ q }|^(ro,T)|{^, q 4)({^, q 4|#t( ro , T ')| {V , 4q}) 
+e - i ^- t '\^ g (r ,r)^l(v ,r')\ip g ) 

x({^ ) q}|*|(ro,r)|{^, q ,})({^, q 4| h(r ,r')m^})] , (22) 

where the sum is over all final states |{</>i', q '}) in the 
excited state manifold. In this expression, we have ne- 
glected contributions involving the product of two cre- 
ation or annihilation operators, a result of the implicit 



rt+At 
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assumption that any atom in the continuum will be re- 
moved from the sample instantaneously. In addition, we 
explicitly carried out the sum over all final states of the 
ground-state field, but not for the excited fields manifold. 
This is because we want to allow for the possibility of se- 
lective detection of the ionized atoms. Following Ref. 
this can be easily achieved by replacing the sum over final 
states in Eq. (|22| ) by a weighted sum 

q'), ( 23 ) 

where 72.(i',q') is the detector sensitivity to atoms in 
state \4>i'^')- In practice, we have in mind energy- 
selective detectors, TZ(i' , q') — *■ TZ(cj), and the degeneracy 
of the levels must then of course be accounted for. 

There is a fundamental distinction between the situ- 
ation at hand and Glauber's photodetection theory, be- 
cause in the present case both the detected and detector 
fields consist of matter waves. There is a complete sym- 
metry between these two fields so far, and their roles are 
interchangeable. In order to break this symmetry and 
to truly construct a detector, we now make a series of 
assumptions on the state of the detector fields ^(r, t). 
Physically, this amounts to making a statement about 
the way the detector is prepared prior to a measurement. 
Specifically, we assume that all atoms are in the ground 
state, vE^ro, 0)|{V'i,q}) = 0, and that any atom in an 
ionized state will be removed from the sample instanta- 
neously, as already mentioned. In that case, the second 
term in Eq. ( p^ ) vanishes. 



the Fourier component of the atomic correlation func- 
tion G^t, i'; r , r ). For stationary fields, the Wiener- 
Khintchine theorem implies that tuning the detector sen- 
sitivity IZ(E) yields the spectrum of the Schrodinger field 
*,(r,0). 

In the case of broadband detection, in contrast, the 
energy distribution AEd of the ionized states is much 
broader than AE g . This situation can be realized e.g. by 
exciting the ground state with a broadband laser pulse 
and detecting the resulting electrons (or ions) with a 
broadband detector TZ(E) ~ constant. Assuming that 
the spectrum of the ground atoms Schrodinger field is 
centered at u> we find 

r bb ~ ri(r )G A (t,t;ro,r ) 7 (25) 

where we have introduced in prevision of the following 
discussion the "detector cross-efficiency" 

7? (r 1) r 2 )=^fi i (r 1 )n*(r 2 ) 




di-(tf j(r 2 , t + r)*|(n, t))e-* a - u ^ T , (26) 



from which the usual detector efficiency is simply recov- 
ered as 77(r ) = 77(1*0, ro). As expected, a broadband 
detector is not able to resolve any spectral feature of 
the Schrodinger field, and only measures the local atomic 
density, like off-resonant imaging. 



2. Energy-selective detectors 



3. Higher-order correlations 



To illustrate this result, we consider the situation of 
energy-selective detectors, and discuss the limits of nar- 
rowband and broadband detection |],||. In the first case 
the detector bandwidth AEd is assumed to be much 
narrower than the energy width AE g of the ground 
state Schrodinger field, which is determined solely by the 
spread in center-of-mass momentum (temperature) since 
all atoms occupy the same internal state. The reverse is 
true in the second case. 

For a narrowband detector, the substitution of Eq. (J22 
into Eq. (p2]) yields readily 



Ai 



r nb (uj)oc / dTe- t{u - UL)T G A {t,t + 



T;r ,r J 



(24) 



Here Tilo is the energy of the registered photoelec- 
trons, and we introduced the ionization rate r nb (to) — 
w n b(to)/At and the normally ordered first-order corre- 
lation function of the ground state Schrodinger field 
G A (t, f ; r , r ) = (*t ( rQj t)* ff (r 0l If)). The only time in- 
tervals that significantly contribute to the integral (|24|) 
are such that At ~ 1/Au) g . For large enough de- 
tection times this integral can therefore safely be ex- 
tended to infinity. In that limit, the detector measures 



The detection of higher-order correlations of the 
Schrodinger field can be achieved by a straightforward 
generalization of the ionization detector. For instance, 
second-order coherence measurements can be carried out 
by focussing the laser at two locations ri and r 2 , in which 
case 

V = hJ2 ^(r^r^^e-^ + H.c. 

/*=1,2 3 

The joint probability to ionize an atom at ri and another 
one at r 2 is then 

/t+At pt+At pt+At pt+At 

dri / dr 2 / dr 3 / dr 4 

" Jt Jt Jt 

xe -iu L{Tl+T2 -r s -r.) n * ( ri )0* 2 (r 2 )%(r 2 )%( ri ) 

x (ri,Ti)^- 2 (r 2 , r 2 )*] 3 (r 2 ,T 3 )*t 4 (r 1) T4)) 

x (*t( ri)Tl )$t (r2)T2 )^( r2) r3 )* g ( ri ,T4)). (27) 

It involves two detected atoms, hence it is now neces- 
sary to properly account for the quantum statistics of 
the measured particles. For this purpose, we describe 



4 



the ionized atoms as ion-electron pairs, whereby the elec- 
trons are described by the creation and annihilation oper- 
ators c ks and Cks satisfying Fermi commutation relations 

[ c ka>cj c v]+ = <5ss"5kk'- Here s labels the electron spin 
and k its momentum. We similarly introduce ion cre- 
ation and annihilation operators aj^aks, also satisfying 
Fermi commutation relations (for bosonic atoms.) For a 
spin-zero atom, the atomic mode operators bj t q can be 
expressed in terms of the ion and electron operators as 



^ q = |jq)<0| = £ |kk' SS ')<kk' SS '|jq)<0| 



kk' 



E 

ks 



(28) 



where cpj (k) are electron wave functions in k-space, Cj s = 
^ k tfj (k)cks and we have assumed that the center-of- 
mass wave function is e lq ' r with r being the ion (or atomic 
center-of-mass) position. Due to spin conservation the 
values of electron and ion spins are clearly opposite. 

Substituting this result into Eq.(p7j) yields, in the case 
of broadband detection, 

/t+At rt+At 
dn J dr 2 

x($t(r 1 ,T 1 )$t(r 2 ,T 2 )$ g (r 2 ,T 2 )$ g (r 1 ,r 1 )) 



dn J dv 



r t+At 



X 



g 

t+At 

+%(ri,r 2 ) / dri 



x 



;*l,( r l' T l) $ U r 2' T l) $ 9( r 2^l)* g (ri,Ti))}, 



(29) 



where the detector sensitivity ^(ri, r 2 ) to processes in- 
volving electron exchange is 



pt+At f*t-\-L\t />t+zir 

r/ x (ri,r 2 )= / dr 2 / dr 3 / dr 4 
Jt Jt Jt 

affKq 

J(w K (Tx— T 3 )+W q (T2— T4)+a>„(Tl — T4,)+Wp(T2 — T 3 )) 

xi^^)! 2 !^^)! 2 ^^)^^)^^)^^) 

+^(r 1 )O Q (r 2 )^(r 2 )fi /3 (r 1 )|^(r 1 )| 2 |0 q (r 2 )| 2 



rt+At 



r t+At 



x e 



"0 

i(tt) re (n— T4,)+W q (T2— T3)+0) a (n— t 3 )+w^(t2— t 4 )) 



(30) 



The first term in Eq. (|29j) is familiar from double pho- 
todetection, with the usual exchange contributions from 
the detected field. The second term is an additional ex- 
change term due to the fact that the detector field is a 
single Schrodinger field. Its origin is the interference of 
the detector field at points ri and r 2 . It is absent in con- 
ventional photodetection theory, a result of the implicit 



assumption that the two detectors used are distinguish- 
able. [J Finally, the term proportional to r\ x results from 
the fact that electrons do not know from which atom 
they originate. We note that these last two terms can 
be eliminated by using a gated detection scheme H that 
eliminates the contribution of the exchange terms in the 
detector field. In practice, such gating can be achieved 
by using non-overlapping short laser pulses to ionize the 
atoms. In that case, the ionization scheme simply yields 
normally ordered correlation functions of the Schrodinger 
field 

G( n \xi, x n ) = (¥( Xl ) . . . &(x n )*(x n ) . . . $(xi)), 

(31) 

in complete analogy with the optical situation. This also 
justifies defining a Schrodinger field as field coherent to 
order Af if its normally ordered correlation functions 
G(") factorize for all n<Af. 



III. EXAMPLE: BOSE-EINSTEIN 
CONDENSATION 

A. Hartree description 

To illustrate the ideas developed in the preceding 
section, we consider a quantum-degenerate iV-particle 
bosonic system described by the state 



i 



N 



d{r i }/ JV ({r J },<)n^t (ri) |o) ; (32) 



where the TV-body wave function /7v({i"i},i) is totally 
symmetric in its arguments. If the sample forms a con- 
densate, it is described to an excellent degree of approxi- 
mation by a Hartree wave function, whereby the iV-body 
wave function f^({ri},t) factorizes as a product of the 
form /Ar({ri},i) = ni</>7v(i"i, t). That is, all atoms in 
the condensate are described by the same Hartree wave 
function (?W(r, t) and the iV-particle state reduces to 



W))n = 



1 



dr^(r,t)*t( r ) | ). (33) 



N 



The equation of motion for ^jv( r ) t) is obtained from the 
variational principle 



1 A similar comment can be made about the position mea- 
surement scheme discussed in Refs. Q|5). In that case, the 
absence of the detector exchange contribution can be traced 
back to the assumption that the set of states excited at 
each location are distinguishable. While this approximation 
is usually reasonable, it becomes questionable in situations 
involving quantum degenerate gases such as Bose-Einstein 
condensates. 
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Mr) 



N 



= o, 



(34) 



where TL is a manybody Hamiltonian. 

We expand then the Schrodinger field as [0 

tft(r) =^(r)ot + 5^t( r ) j 



(35) 



where is the creation operator for a particle with the 
Hartree wave function c/)n(x) and the operator (^T( r ) 
creates particles in all other states. Further assuming 
that the condensate state is stationary 



4> N (r,t) = ( f> N (r)e-^ h 



(36) 



with [in being the chemical potential, the state of the 
system is simply 



-ifiNt/h 



N 



(37) 



where we have used the orthogonality of the condensate 
state to all other modes, i.e. J dr(pN(r)5ij}' (r)|0) = 0. 
Hence the Hartree approximation is equivalent to the as- 
sumption that the condensate is in a number state of the 
self-consistent Hartree "mode." In this description the 
condensate is in a density coherent state (|t]). Indeed, 
\4>(t)) n is easily seen to be an eigenstate of the field den- 
sity operator /5(r), 



p(r)\i J (t)) N = Np(r)m)) N 



(38) 



where p(r) = |</>7v(r)| 2 is the local Hartree density of the 
condensate. It is straightforward to see that this state is 
density-coherent to all orders. However, just like single- 
mode number states of the electromagnetic field, it does 
not exhibit field coherence past first order coherence. 



N N / 



W)) =e-*J2—e- lN " Nt/h ^) N \°)- (39) 



N 



In contrast to a pure condensate state with fixed number 
of atoms in the Hartree ground state, this state can be 
characterized with a nonzero "order parameter" defined 

by | 



$(r,t) = (^(t)|*(r)^(t)) = VN^ R (r)e-^/ h T N (t), 

(40) 

where fx = ^ + N/x^-, v'n = d^N /dN\ N= fj and 



,-2ip'(N-N)t/h 



(41) 



This order parameter exhibits periodical collapses and 
revivals due to the dispersion of the chemical potential 
over the particle number variance. Physically, this is a 
consequence of the fact that while the initial state of the 
condensate is field (Glauber) coherent, it does not remain 
so in the course of time. Hence, the state (|3^) is neither 
field, nor density coherent. 

Interestingly, it turns out that the one-time field coher- 
ence functions measured in ionization experiments (Sec. 
II. C) are factorizable. Indeed, the first order field coher- 
ence function in this case reads 

G«(r;0HiV|^(r)| 2 

and the second order field coherence function is 

G( 2 '(r,r';0) = iV 2 |^(r)| 2 |^(r')| 2 , 

so that the normalized second order coherence function 
ff (2) (0) = L 



B. Wave packet description 

For small systems containing a finite number of atoms 
the assumption that all particles are in the condensate 
state is actually not realistic. Rather, the particle num- 
ber in a condensate fluctuates, so that only the mean 
number of atoms in the condensate is known. In this 
picture it is appropriate to describe the condensate as a 
wave packet in Fock space, and assume that the particle 
number distribution is sharply peaked around N . Fol- 
lowing Ref. Q we assume a Poissonian particle number 
distribution and thus, the state of a condensate reads 



2 Note that this representation of the Schrodinger field op- 
erator is different from the Bogoliubov description, a conse- 
quence of the fact that the system is not taken to be in a field 
coherent state. 



C. Spontaneous symmetry breaking description 

As a final possible description of the condensate, we 
consider the standard spontaneous symmetry breaking 
approach whereby the Schrodinger field operator is re- 
placed by a c- number 

#(r,i) $(r,t). (42) 

This description is equivalent to the assumption that the 
state of the system is an eigenstate of the Schrodinger 
field operator ^f(r,t)\ip) = $(r, t)\if>), or, in other words, 
that the condensate is in a Glauber coherent state. This 
state can be obtained from the previous wave packet de- 
scription in the thermodynamic limit when the dispersion 
of ii over the particle number variance becomes negligible 
and /Ltjv can be approximated by /ijy. In that case 

*(r,t)|V) = y/tie-^^irM (43) 



G 



the condensate remains in a Glauber coherent state at all 
times. It is field coherent, but not density-coherent since 
for that state one finds easily that 

£> (2) (r,r';0) = L>«(r)L>W(r') + £> (1) (r)<5(r - r'). (44) 

These results illustrate how different descriptions 
of the condensate make different implicit assumptions 
about its coherence properties. These models are 
amenable, at least in principle, to experimental tests. We 
note finally that all descriptions reviewed in this section 
become equivalent in the thermodynamic limit as far as 
their coherence properties are concerned, despite the fact 
that their order parameters are different. 



Disregarding mathematical rigor, states which are den- 
sity coherent if all taken at the same time can be 
constructed as follows. Define 



1 



W 



m: • • • j 



Vi) 



Vc 



9Hv m )...&(vi)\0) 



(48) 



with C = (0|#( yi ) . . . % ro )*% m ) . . . *%i)|0)|] This 
state describes a system of m particles which are localized 
at points y\, . . . , y m . For this state the density correlation 
functions are of the form (47) with 



(49) 



IV. MULTIMODE DENSITY CORRELATIONS 
A. Density coherent states 

In this section we further develop the notion of den- 
sity coherence and introduce density-coherent states for 
a multimode Schrodinger field.^j 

^From the definition (]l5|) of the n-th order density cor- 
relation function, we have readily that 



which is a direct consequence of the relation 



D^\x 1 ,...,x n ) = D^(x n ,...,x 1 y 



(45) 



If all taken at the same time D^ n > is real. Similarly 

to the case of field correlation functions Q one can derive 
the inequality 



D ' (xi , . . . , x n , x n , . . . , X\ ) x 
D (2m \y m ,...,y 1 ,y 1 ,...,y m ) > 
\D^ +n \x 1 ,...,x n ,y 1 ,...,y m )\' 



(46) 



The state of a system is said to be Af-th order density 
coherent if all density correlation functions up to order 
Af factorize, i.e., 

D^(x 1 ,...,x n )=D^(x 1 )...D^(x n ), n < TV, (47) 

where we recall that D^(x) is nothing but the expecta- 
tion value of the Schrodinger field density, see Eq. (|ll]) . 



,yi) 



,yi) 



(50) 



Due to dispersion the states \y m , ■ ■ • , J/i) will not remain 
density coherent in the course of timeJj 

For the sake of comparison we also compute the equal- 
time first-order correlation function for a general two- 
particle state 

\4>) = ^= J dx 1 dx2^ux 2 )^(x 1 )^{x 2 )\0) (51) 

with symmetric two-particle wave function 4>(xi 1 X2). 
Denoting the marginal distribution J dx 2 \4i{xi 1 x 2 )\ 2 by 
p(xi) one obtains 

D^( Xl ,x 2 ) = 2\(f>(x 1 ,x 2 )\ 2 + 2p(x 1 )S(x 1 - x 2 ), (52) 

compare with Eq. (p[). This correlation function is not 
factorizable in general. 



B. Thermal fields 

A simple system which can be discussed in the context 
of density correlations is provided by the thermal multi- 
mode Schrodinger field. Here for the sake of simplicity 
we consider an ideal thermal Bose gas of atoms without 
condensate fraction (i.e. above the Bose condensation 



3 The discussion in this section can also be applied to a more 
general class of density correlation functions. To this end, 
consider a one-particle observable B with (discrete or contin- 
uous) eigenvalues b n and eigenvectors \b n ). One might choose 
B — p, for example. In analogy to ^(r, t) and p(r,t) one 
can construct operators 4^(& n ,t) and p(b„,t). In terms of 
these operators one defines the density correlation functions 
according to Eq. ([L5|). The subsequent derivations then apply 
similarly to these correlations functions. In case of B having a 
discrete spectrum the mathematical problems connected with 
the normalizability of the density-coherent states do not arise. 



4 The normalization constant C is not well defined mathe- 
matically, of course. This is due to the continuous nature of 
the eigenvalue spectrum of the position operator f. If one 
considered density correlation functions for a discrete opera- 
tor B (as outlined in the previous footnote) these problems 
would not arise. 

5 However, the momentum eigenstates will be "momentum 
density coherent" at all times for free particles, i.e. if H = 
p 2 /2m. 
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transition temperature T c ) [] in which case the state of 
the system is described by the density operator 

p T = e -^* {ek -^ a * ak/kBT /Z, (53) 

where e k is the eigenenergy of the k— th mode, ks the 
Boltzmann constant, T is temperature, \x the chemical 
potential and Z the partition function. The equal-time 
first-order field coherence function of this system is 

G( 1 \x)=Tr[p T &(x)i>(x)}=J2 ,f k\Mx)\ 2 , (54) 

k 

where we have expanded the Schrodinger field operator 
in terms of mode annihilation operators in the usual way 
as ty(x) = (j>k(x)a' k , and the mean number of atoms 
in mode k is n k — Tr{pTa\a k ). Similarly, the equal-time 
second-order field coherence function is 

G ( - 2) (x 1 ,x 2 ) = Tr[p T &(x 1 )¥(x 2 )i>(x2)i>(x 1 )} 
= ^n kl n k2 (\(j) kl {xi)\ 2 \(t) k2 {x2)\ 2 

+4>* kl (xifyk! {X2)(t>t 2 (* 2 )</»k 2 (xi)) (55) 

These correlation functions are not factorizable. The one- 
point normalized second-order coherence function is 

g [2) (x) = G {2 \x,x)/G {1 \x)G [1) (x) = 2. (56) 

The equal-time second-order density coherence function 
is 

D^(xi,x 2 ) = Tr[p T ^Hxi)^(x 1 )^(x2)^(x 2 ) 

= G^(x 1 ,x 2 )+G^(x 1 )S(x 1 - x 2 ) (57) 

and is likewise not factorizable. Thus the thermal 
Schrodinger field is neither field nor density coherent. 



V. FIELD VERSUS DENSITY COHERENCE IN A 
DYNAMICAL SYSTEM: THE 
BINARY-COLLISIONS ATOM LASER 

In order to illustrate how differently the field and den- 
sity coherence of a Schrodinger field can behave in a dy- 
namical system, we compute these properties for a model 
of a binary collisions atom laser. We show in particular 
that elastic collisions, while being extremely detrimen- 
tal to the field coherence, and hence the linewidth of the 
laser, have almost no influence on its density coherence. 



KiN t 










K 2 


Ki(N+l) 




%211 




K 



FIG. 1. Schematic three-mode atom laser model. 

This binary collisions atom laser has been investigated 
in various publications |ll|-[lrj] . One considers a resonator 
for atoms, realized e.g. by optical fields, in which only 
three atomic center-of-mass modes are taken into account 
explicitly (cf. Fig. [I]). Bosonic atoms in their ground elec- 
tronic state are pumped into an atomic resonator level of 
"intermediary" energy (mode 1). They then undergo bi- 
nary collisions which take one of the atoms involved to 
the tightly bound laser mode 0, whereas the other one 
is transferred to the heavily damped loss mode 2. This 
latter atom leaves the resonator quickly, thereby provid- 
ing the irreversibility of the pumping process. A macro- 
scopic population of the laser mode can build up as soon 
as the influx of atoms due to pumping compensates for 
the losses induced by the damping. 

For the master equation of this atom laser model one 
makes the ansatz 

W = -i[H + B C ou W] + K a V[a a ]W + k x (JV + l)V[a x ]W 
+K l NV[a\]W + K 2 V[a 2 ]W (58) 

with h = 1. In this equation, the second quantized for- 
malism is used in which each center-of-mass atomic mode 
is associated with an annihilation operator a,, and W 
denotes the atomic density operator^. The free Hamilto- 
nian is given by 

H o = ^ Wiffltoi, 

i=0,l,2 

uii being the mode frequencies. For the following we set 
luq + uj 2 = 2uj\. This condition can be fulfilled, e.g., if 
the atomic cavity is realized by time-modulated optical 
fields ' 

The operation of the atom laser relies on binary col- 
lisions between the atoms in the resonator. The general 
form of the corresponding interaction Hamiltonian is 



6 The calculation of the second-order field coherence function 
for the thermal Bose gas with both interactions and conden- 
sate fraction present can be found in [10| . 



7 Since we consider ground state atoms only, they are fully 
described by their center-of-mass quantum numbers. 
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Hr, 



i<j,k<l 



V ijk ia\a\a k ai 



(59) 



where Vijki are the matrix elements of the two-body in- 
teraction responsible for the collisions. However, for the 
present investigation we restrict our attention to the sim- 
plified form 

H co i = V Q2 iial a 2 a i a ^ + ^1102 <A a\ ao a>2 + Vboooao a o a o«o 

(60) 

in which (besides the pumping collisions) only collisions 
involving ground state atoms are retained. The damping 
rates of the cavity modes are given by the coefficients 
Ki, and the strength of the external pumping of mode 1 
is characterized by the parameter N, which is the mean 
number of atoms to which mode 1 would equilibrate in 
the absence of collisions. The superoperator T> is of the 
Lindblad form and is defined by 



V[c]P = cPc f - \{ScP + Pch) 



(61) 



with arbitrary operators c and P. 

In order to achieve a sufficiently high degree of irre- 
versibility it is necessary that ki is much larger than the 
damping rates of the other modes. This suggests to adi- 
abatically eliminate this mode, an approximation that 
leads to the simplified master equation JTll— 



(62) 



-i[H c ,p] + n V[a ]p + Kl {N + l)V[ ai ] P 
-KiNT>[a\]p + TV[ala% 



Equation (|6^) is written in the interaction picture with 
respect to Hq = uJoaQaQ+ujia\ai, and the reduced density 
matrix p is p = Tr mQc j e 2 [IV]- The reduced collision 
Hamiltonian H c is 



H c — Vboooao a o a o a o, 



(63) 



and T = 4|Vo2ii| 2 /k. Consistently with Ref. |T^] we call 
the limiting cases T <ti kq and r > (to the weak and 
strong pumping regimes, respectively. 



A. Linearized fluctuation analysis 

In order to obtain analytical approximations for the 
correlation functions 



and 



Gf ) (r) = (at(0)a].(r)a i (r)a i (0)) 



Df\r) = (a](r) aj (T)a}(0)a 3 (0)) 



(64) 



(65) 



in the two-mode sytem a linearized fluctuation analysis 
can be performed PJl^j. To this end the master equa- 
tion (62) is converted to a Fokker-Planck equation us- 
ing the P-function representation as described in p2[. 



This equation can be transformed to polar coordinates 
a j = V™j elc ^ j > where- oi j denotes the complex ampli- 
tudes originally appearing in the Fokker-Planck equation 
E^Eg) . This leads to stochastic differential equations 



dn = [Tnl(n + 1) - n n ]dt + dS no , 
dm = [«;i(iV - ni) - 2Tn\(n Q + l)]dt 
d(p = [-Vbooo(2n - 1) - Vbioini]^ - 
d(j) t = [-Viiu (2ni - 1) - V 010 in ]dt - 

The correlation matrix D for the stochastic forces dS T 
(dS no , dS ni , dS<j, , dSfa ) is given by 



I- dSr, 
dSfc 
dSs, 



(66) 
(67) 
(68) 
(69) 



/ 2Tnln 
-2Tn\n 
- 2Voooo 7 io 
\ -Vbioi^o 



-2Tn\n 
2mNni ~ 2Tn\n 
-Voicn^i 
-2Vimni 



— 2Vbooo^o 
-Voioim 
r7i2/(2n ) 

rm/2 



-Vbioi«o \ 
-2Viiii?ii 
rm/2 

KiN , Trig 
2ni 2 / 

(70) 



In the limit n ^> 1 one obtains from Eqs. ( pq ) and 
( |67| ) the above-threshold semiclassical steady-state pop- 
ulations |l2| 



m = 




(71) 
(72) 



the threshold condition being N > ^/ko/T. The drift 
terms in Eqs. ( |66| ) - ( |69| ) and the correlation matrix D do 
not depend on the phases <pj. This means that the time 
evolution of the atom numbers rij is not influenced by 
the phase dynamics and is thus completely determined by 
Eqs. ( |66|) and (|67j ) alone. To proceed further we introduce 
the fluctuation variables Srij = rij — rij. In the linear 
approximation their time evolution is given by 



dSn 



-kSndt + dS r , 



(73) 



where Sn T — (Sno, 5ni) and the matrix k is obtained by 
linearizing the drift terms in Eqs. (^6|) and ( |67j ) around 
the steady-state values rij. The correlation matrix D„ 
for the stochastic forces dS ^ = (dS na , dS ni ) is given by 
the upper left 2 x 2-minor of the matrix D of Eq. (170) 
after replacing rij by hj . 



B. Second-order correlation function G' 2 ' 

In the linear approximation the steady-state second- 
order correlation function for mode j is given by [pLR 



(74) 



(t > 0) where <x is the steady-state covariance matrix 
AD„ + (k-£)D„(k-£) T 



2SA 



(75) 



4UU 



o 
o 






'in 



44 



4k + Kini/no ' 
qi = «i + 4K «o/ni 



(80) 
(81) 



Their inverses can be interpreted as the relevant 
timescales for the dynamics of the atom number fluctua- 
tions. Numerical comparison to Eq. ( ff4|) shows that the 
approximation (|7^) is very accurate, in general. The time 

dependence of the second-order correlation functions is 

(2) 

thus purely exponential. The correlation function G 

decays on a timescale of the order of kq whereas the time 
(2) 

evolution of G\ is much more rapid. 



C. Density-correlation function 



FIG. 2. Dependence of (Too on no- The curves are cal- 
culated using Eq. ( fF^ ) for parameter values F = 0.07ko (as- 
cending curves) and T = 15ko (descending curves), n\ = 20reo 
(full), ki = lOO^o (dashed). The single points show numeri- 
cal results for F = 0.07ko, k i = 20fto (•) and Ki = IOOko (*), 
respectively. 



Using Eq. (10.5.28) of Ref. |§ we find that in the 

linear approximation the steady-state density-correlation 

(2) 

function D) for mode j is given by 



Df\T)=n J (e-n n +GY>(T). 



-k-n 



(82) 



with A = det k and E = Tr k (we use the numbers 0,1 as 
indices for the 2 x 2-matrices). For the matrix elements 
of <j one obtains the expressions 



coo 



[n„(2Ko - ko/«i) + n (Kini + k + ki) 
+K 2 1 n 1 / (4K )]/(4rnin + ti), 
Ko(2ni — V)fiQ + Kin\ + Korii 



4k "-o/^i 
f 01 = cr w = -ni/2. 



(76) 
(77) 
(78) 



In the far-above-threshold limit uq 1 the covari- 
ances can be approximated as er o ~ 3(^1 — ^)no an d 
o"xi ~ gC^i - 5)^1) respectively. For ni > 1/2, i. e. in the 
weak-pumping regime, the second-order correlation func- 
tions thus show bunching, whereas for fi\ < 1/2, i. e. in 
the strong-pumping regime, they are anti-bunched. The 



normalized second-order correlation (^(O) = G ( o z '(0)/nQ 

(2) 

goes to 1 with 1/rio for large uq. In contrast, g\ '(0) — > 
3/2 — l/(4ni). For the case of fix 3> 1 these results are 
in agreement with the conclusions of Ref. Eh To illus- 
trate the physical contents of Eqs. (|7^) and Q77[) in Fig. || 
the dependence of <Too on uq is shown for the weak and 
strong-pumping regime and different values of K\ . 

An expansion of e~ kT in the parameter n^/ni yields 
the approximate result 



,(2) 



Gf\r) 



-qjT 



(79) 



for the time dependence of the correlation functions (74) 
Thereby, the qj are the two eigenvalues of the matrix k 
For fio/ni > 1 they are approximately given by 



From Eqs. @ and @ it follows that Df ] (0) > n|. In 
analogy to Eq. ( |7^ ) one obtains the expression 



- qj i 



(83) 



(2) 

for the explicit time dependence of -D) . The inten- 
sity fluctuation spectrum, i.e. the Fourier transform of 
Dj(r) — fij, is thus a Lorentzian. 

It should be noted that all results of this section are 
independent of the rate of elastic collisions (as quanti- 
fied by Voooo) between atoms in the laser mode. This 
means that in the two-mode description the atom laser is 

(2) 

second-order coherent (at least in the sense of <?q (t)»1 
for all t) even if Voooo is large. In such cases the first- 
order correlation function G' 1 ' will decay very rapidly so 
that the laser is not first-order coherent, and is charac- 
terized by a large linewidth. 



D. Numerical examples 



As discussed in detail in Ref. |15| the master equations 
(f58| ) and ( |62|) may be solved numerically with the help 
of quantum Monte Carlo techniques |^,|l| . The correla- 
tion functions G^ and D^> can be calculated within this 
approach according to the description given in Ref. pl| . 
In the following we compare some numerical results with 
the analytical predictions of Sees. IV. B and C. 

Strong-pumping regime. Figure p|(a) depicts a typi- 

(2) (2) 

cal result for the calculation of D\ and G\ in the 
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FIG. 3. Numerical calculation of D<, 2) (r) and (t) for 
parameter values F = 15ko, Ki = 20«o, iV = 1.2 (a) and 
r = 0.07so, Ki = 20ko, AT = 3.9 (b). Dashed curves: expo- 
nential decay with corresponding time constants 1/qo- 

strong-pumping regime. In this example, the parameters 
r = 15ko, ki = 20ko, and N = 1.2 were chosen yielding 
a numerical equilibrium population of no = 9.0 (9.4 ana- 
lytically). As should be expected, D ( Q 2) (0)-G ( 2) {0) ps n . 

(2) 

The behavior of L» ' (t) is very well approximated by an 
exponential decay with time constant 1/qo, as predicted 
in Eq. (^||) . From the behavior of G Q 2 ^ (r) it can be in- 

ferred that (t) — l| <C l for all r. However, even 
after a very large number of Monte Carlo simulations 
numerical noise prevents any further details of the time 

(2) 

dependence of G (r) to be identified. This observa- 
tion applies to most calculations in the strong-pumping 
regime. In particular, it could not be unambiguously de- 
termined whether antibunching actually occurs for larger 
values of no- On the other hand, these numerical results 
are compatible with the fact that Eq. ( |76| ) predicts a 
small value of |<?o (r) — l| (compare with Fig. ||). 

It should be noted that a small amount of antibunching 

(2) 

in G Q could be observed in calculations for the three- 
mode model in cases in which Vbooo is large compared to 
Vo2ii- Under these conditions, however, the main effect 
to be observed is a significant decrease in the equilibrium 
population Hq. 

Weak-pumping regime. In contrast to the previous 



t(2). 



11.0 - 




FIG. 4. Numerical calculation of G 1 ' 2 (t) for parameter val- 
ues F = 0.07«o, ki = 20«o, N = 3.9. Dashed curve: expo- 
nential decay with time constants 1/qi- 

case, for the weak-pumping regime it can be expected on 

(2) 

the grounds of the linear analysis that g (0) is signifi- 
cantly different from 1 if no is not too large (cf. Fig. |^). 
As exemplified in Fig. |^(b) this prediction is indeed con- 
firmed by the numerical calculations. There the correla- 

(2) ^(2) 

tion functions D Q and G are shown for the parameter 
values r = 0.07koj k i = 20ko> and N — 3.9. Unfortu- 
nately, the quantitative results of the linear analysis in 
the weak-pumping regime are not very accurate for low 
values of no (to which the numerical computations have 
to be restricted due to time constraints). For example, 
for the above parameters Eq. ( |7l| ) yields no = 1.2 which is 
much smaller than the numerical value of 11. However, 
Eq. (|8(]) still constitutes a good approximation to the 
decay rate of the correlation functions if it is used with 
the numerically determined values of no and fi\. This 
is demonstrated by the dashed curve in Fig. ||(b) which 
depicts an exponential decay with a time constant calcu- 
lated in this way. A similar agreement was also found in 
other examples. 

In Fig. H the results of several numerical calculations 
of <Too are shown as a function of the numerical value of 
no for the parameters T — 0.07ko, k\ = 20ko (•) and 
Ki = 100ko (o). The numerical values of er o should be 
understood as having an error margin of at least ±10 - 
15%. The results depicted agree in order of magnitude 
with the predictions of Eq. (76) and also demonstrate a 
dependence of <r o on uq and «i similar to the analytical 
one. 

Second- order correlation function for the pumping 
mode. An example of the behavior of the second-order 

(2) 

correlation function G\ (t) in the weak-pumping limit is 
shown in Fig. |]. There, the parameter values T = 0.07koj 
Ki = 20kq, and N = 3.9 were used. The order of mag- 

(2) 

nitude of G\ '(0) as well as the temporal decay rate of 

(2) (2) 

G\ (which is much larger than the decay rate of G ) 
are in good agreement with the analytical predictions. 
It should be noted that the time evolution of the mean 
population of the pumping mode (a\ai(t)) starting from 
an initial vacuum state contains both characteristic time 
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scales and 1/qi- This is quite different from the be- 

(2) 

havior of G\ which is characterized by 1/qi alone. In 

(2) 

the strong pumping limit the time dependence of G\ 
could again not be recognized due to numerical noise. 
This is consistent with the fact that the linearized fluc- 

(2) 

tuation analysis predicts a small value of \g\ (0) — 1| for 
the parameter values investigated. Furthermore, the cal- 
culations indicated that Eq. (|7|) ceases to be valid for 
large T. It was not possible to reduce h\ to very small 
values in which case a more pronounced anti-bunching 
would be expected. 

In conclusion, we see that the results obtained from 
the linear fluctuation analysis describe well the essential 
aspects of the behavior of the correlation functions Dq 

(2) 

and Gq and may be used as a first quantitative estimate. 



VI. CONCLUSION AND OUTLOOK 



In this paper, we have adopted an operational ap- 
proach to introduce several classes of coherence of the 
Schrodinger field. Of particular importance are den- 
sity coherence, which is connected to far-off resonance 
imaging measurements, and field coherence, which can 
be measured in ionization-type measurements. One can 
readily imagine further classes of coherence associated 
with other types of measurements, but they are proba- 
bly not as important as field and density coherence. 

One question of considerable significance in the future 
will be to quantify the usefulness of various sources of 
Schrodinger fields for specific applications. In optics, 
one of the most important characteristic of lasers is their 
spectral width, and higher-order coherence plays a lim- 
ited role in most cases. By analogy, past studies of atom 
lasers have concentrated on their spectral width, as de- 
termined by their first-order field correlations. It as been 
found that this linewidth can be quite broad, specially 
in the presence of elastic collisions. Indeed, things can 
be so bad that the atom laser linewidth is broader than 
the natural linewidth of the atom cavity, in sharp con- 
trast to the optical situation where, of course, the reverse 
is true. It is not clear however whether this is a useful 
way to determine the quality and usefulness of an atom 
laser. While this is likely to be the case in some inter- 
fcromctric applications, other possible uses of atom laser 
beams, such as coherent lithography, may well require 
only a high degree of density coherence, in which 
we have seen, elastic collisions do not play a detrimental 
role. Hence, we believe that it is important at this point 
to start analyzing in detail the coherence requirements of 
specific atom laser applications, so as to optimize their 
design. 
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